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flj ', Abstract. We provide a concrete and systematic connection between the statistical 

Ch ' physics of the Ising ferromagnet on a Cayley tree, which is an old solved problem, and 

the study of memory in exponentially expanding spaces, which is of current interest. 
C^ ' We provide the order parameters which clearly signal the (known) phase transitions. 

C/3 ' The correlators are found in a new simple way, which may be of interest for other tree 

like systems, and which allows the computation of higher order expectation values of 

ethe order parameter. Numerical Monte Carlo simulations are used to make transparent 
the existence of the two different transition temperatures. 
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1. Introduction 

The objective of tliis paper is to make a concrete and transparent connection between 
two separate fields of study that are interested in the same problem [H |2l |3l HI |5l |6] and 
to bring to the fore old work, which is of very current interest, in modern language from 
a contemporary perspective. 

Our motivation to revisit an old statistical physics problem came from recent works 
showing the natural appearance of tree like structures in certain gravity phenomena, like 
bubble nucleation in eternal inflation [7] or in the near horizon region of black holes [8] . 
The Ising model in the Cayley tree is complicated enough that some interesting questions 
that arise here can be meaningfully posed, but still simple enough to answer all of them 
in a clear way. 

In particular, in this article we develop the connection between questions related to 
memory, which have their root in the branching Markov process approach, and which 
were raised in a different framework by Ref.[9], and the more usual static approach, 
dealing with partition functions and correlators. Building on old knowledge of this 
system and with the aid of new techniques to compute certain correlation functions in 
the Cayley tree, we obtain a clear and simple way the two known critical temperatures, 
the so called Bethe-Peierels [1] and spin glass [10] critical temperatures. Furthermore, 
we point out that the expectation value of the product of the spin at the root and 
the total magnetization at the boundary {a^Mn) is the appropriate order parameter to 
signal Tbp, and ^"^"^ "^ where Vt = ^J {{aQMnY) — (aoMn)"^ is the RMS deviation, is the 
order parameter for Tsg - facts that appear to be not widely known. We also provide 
a rigorous derivation of the Markov transition matrix from the static approach. Also, 
we point out that this conceptual framework can be easily extended to the spin-glass 
and Heisenb erg-type variants of the Hamiltonian - a discussion of some consequences is 
included. 

Monte Carlo simulations on these type of systems in which the boundary is of the 
same size as the system itself have always been problematic, due to subtleties with the 
thermodynamic limit. With the help of this new analytic understanding we are able to 
perform Monte Carlo simulations, showing clearly the two transition temperatures. 

We make clear right at the outset that we consider here the Cayley tree and not 
the Bethe lattice, the difference between which two has been a source of confusion in 
the literature - look for example at Ostilli [TT] . 

2. The Static Approach: Partition Function 

In 1974, Eggarter [1] provided a beautiful way of computing the partition function of 
the Ising model on the Cayley tree. The method is based on two important properties. 
The first one is that, generically, the Ising model Hamiltonian is just a sum of 'bonds', 
i.e, a sum of products of adjacent spins. The second is that the Cayley tree has no loops 
and, therefore, the bonds provide an equally good basis in which to define the partition 



Memory as order-parameter on a Cayley tree 




Figure 1. Figure shows a Cayley tree of depth n ~ i and coordination number z = 3. 



function. The advantage is that the Ising Hamiltonian decouples in these variables. 
More concretely, we start with the following Hamiltonian: 

"H = -J^CTiO-j, (1) 

where the sum involves pairs of spins which are adjacent on the tree. Instead of the 
spin variables, we choose the spin at the root of the tree (Tq, and the bonds 6ij = aiaj to 
define the configuration of our system. The 9 variables can take values ±1, which make 
them effectively spin variables too. The Hamiltonian then takes the following simple 
form 

n = -jY.9,,. (2) 

We see that the Ising model on the Cayley tree is essentially a sum of independent spins 
(the bonds) subjected to a constant magnetic field; in the computation of the partition 
function, the only difference with this dual system is a factor of two coming from the 
sum over the possible states of the spin at the root. Therefore we arrive at [1]: 

Z{J)cayley = 2(2 COsh(/3 J))^^ = 2(2 COsh(/3 J))^"\ (3) 

where A''^ is the number of bonds and A^ is the number of spins. 

With this formalism, namely the 9 variables, the computation of the correlators is 
straightforward. Firstly, odd correlators vanish due to the spin symmetry of the model. 
For the even case we begin with the two point function. We would like to write the spin 
product in terms of 9 variables. To achieve this, we just need to find the path between 
the two spins. This path is unique because there are no closed paths in a Cayley tree. 
Imagine that the path goes over era, ctj,, ■ ■ ■ , ac. We can then write the product as 



(TiCT 



i^] 



aiaaaaabO-b 



O'cO'cO'j 



9ia9a 



■9 



CJ- 



(4) 
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The important point now is that, as commented before, the 9 variables are not coupled 
to each other, and are thus independent variables obeying the same Hamiltonian. 
Therefore, 

{a,a,) = {6, J,, ■■■0cj) = {0.a) (Oab) ■ ■ ■ (Oc,) = (^)"- , (5) 

where dij is the discrete natural distance on the tree, and ^ is a spin variable obeying 
the following simple Hamiltonian 

He = -Je (6) 

with partition function 

ze = 2cosh(/3J), (7) 

and whose one-point function is given by 

<"> = ^^ = ^^r^ = t-h(/JJ). (8) 

Therefore we have, with the defintion a = tanh(/3J) which will prove useful later, 

{(x.aj) = tanh(/3J)'^'^ = a'^^^ . (9) 

Going to higher order correlation functions is also possible. The above technique 
allows us to obtain the general result for a 2ra point function, namely: 

{ai.ai, ■ ■■ai,^) = (tanh(/3J))[^(^s) '^«a>s1™ (10) 

where di^ig are the distances between pairs and the minimization is carried out over 
all possible ways in which n-pairs can be formed from 2n points. A simple geometric 
construction to find this minimum is to select n-pairs at random, join all the paths 
between pairs, and delink any portions on the tree that are traversed an even number 
of times. The newly formed pairs will give a minimum for the sum. 



3. The Dynamic Approach: Evolution of a Markov Chain 

Now we consider the type of approach adopted in Refs. [6l |9]. The problem is the 
following. Consider an Ising spin at time t = 0, which can be up or down with 
probabilites Pq = (pJ,Pq). At time ti this spin gives rise to z spins, each of which then 
gives rise to {z — 1) spins at the next time instant t2, each of which in turn gives rise 
to (^ — 1) spins, and so on. There is an obvious tree structure appearing here, with the 
time t playing the role of the lattice depth in the past static scenario. Now the dynamics 
is determined by a Markov rule: the probability distribution of a spin at level t„+i is 
completely characterized by the one of its parent spin at level t„. This process defines a 
Markov chain, in particular a Branching Markov chain, due to the exponential growth 
of the number of spins at every time step. A Markov chain is defined by its transition 
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matrix. In our ca se, taking into account that the probabihty ratios should correspond 
to the ones coming from the Gibbs distribution we have the following transition matrix: 

^("'"') = 2cosh(/3J)=^r ^''^ 

The Markov dynamics is then defined as 

Pn+m = G Pn- (12) 

A rigorous technique, which uses the above 9 variables, to obtain the transition 
matrix, is the following ( for clarity we define the transition matrix between the spin at 
the root and one of its descendants but the computation can be done similarly for any 
pair of spins) : 

The joint probability density of the spin at the root and one of its adjacent spins 
is given by 

^_ P^PH pjpeoi 
P{cro,cr,)= Yl -^ = -^r^^ee-'r^-' (13) 

pJl39oi pJPdoi 



Z "^ 2ze 

and the reduced probability for the spin at the root is just (as expected) 

The transition matrix is, by definition, the conditional probability P(cri|cro), and so 

G (71, CTo = p. X = = , (15) 

P(o-o) ^e ze 

which is the transition matrix of EqJTTl Now, in order to compute the correlations, we 

first diagonalize the matrix G. It is easy to check that, with a = tanh(/3J), 

(16) 




(17) 

and as expected, no matter the initial distribution, the probability distribution at late 
times tends to the stationary distribution P = (1/2, 1/2). 
The correlator between any two spins is given by 

(or-cTj.) = ^ aiajP{ai,aj) (18) 

o'i,(Tj=:izl 

where P{ai,(jj) is the probability of finding ai,(jj. This probability can be found by 
going to the first common ancestor of the two spins, 0"r-. It is given by 

P{a,,a,)= Y^ Pi(T,\ar)P{a^\ar)P.^ (19) 

(Tr=ltl 



Memory as order-parameter on a Cayley tree 





(b) 

Figure 2. Monte Carlo simulation of Cayley trees with coordination number z = 3. 
Figure (a) obtains the transition temperature popularly called Bethe-Peirels Tbp- 
Figure (b) obtains the so-called spin-glass transition temperature [10]; Tsg deeper 
down in the ordered phase. [Appendix C has some details of the simulation |Appendix D| 
shows why the second-transition temperature carries the name 'spin-glass': quite 
remarkably, the same transition temperature appears when one considers the ± J spin 
glass on a Cayley tree. 



Taking into account that individual spins, and in particular the common ancestor, have 
stationary distributions in the thermal ensemble, and that the conditional probabilities 
P{(7i\ar) are just given by the elements of the matrix C^^^ we arrive finally at the desired 
result 



(a,a,) = a"^- (20) 

So the Markov chain, which is a time process, gives the same value for the correlators 
between the spins as the common static procedures on a thermal tree. This observation 
may be interesting for the model studied in Ref.[7], which, in this perspective, may be 
expected to be equal to some statistical model on a thermal tree. Several observations 
done in that reference become very natural in this perspective (like the p-adic structure 
of the boundary correlators). We leave these interesting issues for the near future. 



4. The two phase transitions 

In the Markov chain approach, a very natural question to ask is how correlated some 
random variable at the n^^ time-step is with respect to the root spin; in other words, how 
much memory has that random variable about its earliest ancestor. For example, the 
correlator {ctoCTj) dies for large distances, no matter the value of f3J. So, as commented 
in Ref. [9], there is no memory for these two random variables. We now show that, 
in fact, the consideration of the correlation between the sum of the spins at the n^^ 
time-step with respect to the root-spin, i.e. (aoMn), is more useful and gives the order- 
parameter for the so-called 'Bethe-Peirels' transition. Furthermore, we will show that 
the consideration of the allied quantity ^'^°^' "' where f2 = -^((o"oMn)2) — (o"oM„)2 is 
the RMS deviation, gives the so-called spin glass transition temperature. Appendix A 
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carries a discussion of the motivation behind studying these two quantities from a 
generic perspective, and in particular its conceptual relation with the mutual information 
between the two random variables. 

Indeed, that the Ising ferromagnet on a Cayley tree has two distinct phase 
transitions is long known in the statistical physics community - this information is 
buried in the papers of Eggarter [1], Matsuda [2], and particularly, the papers of J von 
Heimburg and Thomas 0], and Morita and Horiguchi [3]; however, we believe that our 
paper shows the correct order-parameters for these two transitions for the first time. 

Computing the correlation between the root spin and the sum of the spins at the 
^th igygi ^g obtain: 

(cToM„)= Y. K^.) = iVntanh(/3Jr, (21) 

where Nn is equal to the number of spins in the n*^ layer: Nn = z{z — l)'^"^, where z is 
the coordination number for the type of Cayley tree shown in Fig{TJ Therefore, defining 
7 = 2; — 1, we have 

(cToM„) = (7tanh(/3J))"(l + -) (22) 

7 

It is clear that the behaviour of (aoMn) for 'long times' or, equivalently, in the 
'thermodynamic limit' n — )■ 00, changes completely depending on whether 7tanh(/3J) 
is > or < 1. So, 7tanh(/3J) = 1 signals a transition which is understood to be the 
so-called Bethe-Peirels critical temperature Tbp = — t^ttxt- FigJ2]^£'') shows that this 
phase transition can be obtained from Monte Carlo simulations of finite-systems. 

For computing the deviation of this random variable we first need to compute the 
second order moment. It is given by 

= (7 + l)7"-^ + 2^a'^'^ (23) 

The last sum above was done in Ref.[3]. We compute it in a different way. In particular 
we can derive a recursion relation for the sum, the details of which are provided in 
Appendix B[ We collect here the result: 



Ea-~-iV„( 'V;>-";:7 4). (24) 

where x = 70^ = 7tanh (/3J), and Nn = (7 + 1)7"""^ is, once again, the number of spins 
in the n^^ layer. With the help of Eq JC.lt and some algebra, we can now write down 
the mean square deviation: 

n^ = {{a^Mnf) - {aoMn)^ = iV„(l + ^ ~ ^ "^^ ~ "^ - -x") (25) 

7 X — 1 7 

From this expression one can easily check that the ratio (croM„)/r2 defines two 
phases delimited by the second critical temperature, which is given by x = 1 (Tsg = 
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- — j^ j_ ). In the high temperature hmit, x < 1 this ratio dies in the thermodynamic 

hmit, signalhng the absence of memory as discussed previously. In the low temperature 
regime, x > 1, the ratio begins to grow, diverging in the strict zero temperature limit. 
Therefore, this quantity (o"oM„)/r2 provides a good order parameter on the Cayley 
tree. Figl2]^b) obtains this transition temperature by studying the quantity Vt^ scaled 



appropriately. Appendix C has a discussion of why this quantity is more convenient for 



finite-sized systems. The reason the second transition temperature is called 'spin-glass' 



is explained in Appendix D Quite remarkably, the methods of Eggarter and others carry 



right through when one considers a ± J disordered Hamiltonian on the Cayley tree, and 
there occurs only one transition temperature at precisely - — _'(.j_ , thus justifying the 



name 'spin-glass'. However, look at Appendix D for a discussion on how there is not a 



regular type of 'spin-glass' phase below Tsg- 

One could expect that higher moments of cr^Mn would give us new transition 
temperatures, which would presumably be related to the infinite set of critical 
temperatures found for this system in [5] . With the aid of the new recursive techniques 



discussed in Appendix B, we were able to compute analytically the next moment, 
without finding any new critical temperature. We leave this interesting connection 
for future study. 

We note that the techniques of Ref.[T2] can be used in conjunction with the 
above insights to study the classical Heisenberg ferromagnet on a Cayley Tree. For 
this case, the function that we have defined as a = tanh/3J would merely change to 
ttHeisenberg = coth.(3J — jj , and thus the two transition temperatures would be given by 
10, Heisenberg = I5 and Id Heisenberg — 1' closcly aualogous to the Isiug casc. Indeed the 
same logic should hold for arbitrary m-component vectors though a closed form solution 
is unlikely. 

5. Conclusions and Outlook 

We have shown that the dynamic time Markov chain approach gives exactly the same 
results as the static ensemble approach of statistical mechanics, thus linking two different 
communities. We have obtained the two transition temperatures present in this model 
with new recursive techniques, which may be interesting for other tree like models. 

We have highlighted that coMn is a convenient order parameter for the system, 
which can at the same time be thought of as a kind of memory in the time-picture. 
This quantity signals the Bethe-Peirels transition temperature Tbp- Furthermore, we 
have shown that ^"^"^ "' , where Vl = ^J {{aQMnY) — (aoMn)"^ is the RMS deviation, is 
the order parameter for Tsg- The analytical understanding allows us to tune our Monte 
Carlo simulations, which are usually problematic in these kinds of lattices due to the 
difficulties appearing when taking the thermodynamic limit, to study the right quantities 
showing clearly the transition temperatures. 

The analogue to (JoM„ in other expanding structures such Euclidean Anti de Sitter 
or de Sitter space-times is a very natural and known quantity, the study of which 
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in this light we postpone for the future. Indeed, the fact that in the Cayley tree 
simple correlators provide the same physical results as more involved memory type 
questions concerning the mutual information explains the relation signalled in Ref. [9] 
between quantization schemes and the persistence of memory. The rate of falloff of the 
correlators, determined in our case by the temperature and in the EAdS and dS cases 
by the dimension of the operator being examined, plays the central role in this problem. 
It would also be interesting to study the correlations between the central node and 
the sum on the n^^layer in quantum models; the entanglement between the central node 
and the sum on the n^^layer could be another intriguing direction to pursue. 
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Appendix A. Obtaining the order-parameters 

Suppose we want to measure how correlated two random variables A, B are. The 
quantity that is commonly considered in the Markov process community is mutual 
information, which is defined as 

I{A,B) = S{A) + S{B)-S{A,B), (A.l) 

where 

A 

S{B) =-^P{B)\og{P{B)) and 

B 

S{A, B) = -J2 PiA B) log(P(A, B)) (A.2) 

A,B 

are the usual Shannon entropies. This quantity is bounded by zero from below (This 
happens when P{A,B) = P{A)P{B), which implies that the correlator of the product 
factorizes (AB) = {A){B)). From above, it is bounded by the entropy of the random 
variable with lower entropy. This happens when knowing one of them implies knowing 
the other(for this to happen we need that the entropy stored in the first is at least 
as much as the entropy stored in the second). For this case, one would expect that 
the quantity (AB) — {A){B), which is just the often-studied correlation function in 
statistical physics, becomes large. However, the ratio - — '~^ '^ ' , where aAs is the 
standard deviation, must also be studied because we may have a scenario where cab is 
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also large, and thus the knowledge of one of the two random variables A leads to very 
little predictive power for the variable B even though the two have a large correlation 
{AB). These two different types of quantities gives us order-parameters that mark 
precisely the two different phase transitions discussed in the main paper. 

Appendix B. Recursion Relations to compute deviation 

Here we obtain the sum X]i<7'^*^' where all z,j are points on the n^'^-layer boundary, 
with the help of a recursion relation that is obtained by means of a graphical iteration 
which constructs the tree step by step, from the boundary to the root. We begin by 
considering two spins at the boundary at distance 2 in tree units (the minimal distance). 
At this initial step the sum is, since there are (2) ways of choosing the two spins, just 
trivially (2)0^ and which provides an initial condition for the iteration. Now we go one 
step higher, by stitching together 7 of these 1-layer trees to one spin above, and thus 
also connecting spins which are at distance 4. The sum in this case is found to be 7 
times the first sum, corresponding to the 7 cases in which the two spins lie within the 
same old branch, and (2)7^ terms with value equal to a^, which give us the cases in 
which one spin lies in one old branch and the second spin in another. This pattern 
is easily seen to hold until one step before the last iteration (the last iteration should 
include one branch more in the tree we are studying, and has to be treated separately). 
Calling this iteration Sk the recursion relation takes the following form: 

5.+i = 75.+ (])fV(^^^\ (B.l) 

where 5*1 = (2)^^^- This can be seen to be a geometric sum by the following 
transformation Sk — ;• (2)7^^^C*a;- The solution is then 

S, = 1'-'- r^, (B.2) 

X — 1 

where x = ^0? . 

Adding 7 + 1 braches to form the top root-node (because of the type of Cayley tree 

we have considered), we have 

^-^ 27 X — 1 2 

where A„ = (7 + 1)7""^ is the number of spins in the last layer. 

Appendix C. Monte Carlo Simulation 

Here we show how to exploit the analytical understanding to obtain the transition 
temperatures in the thermodynamic limit, with the help of Monte-Carlo simulation on 
finite-sized systems. We showed in the main paper that 

(aoM„) = (7tanh(/3J))"(l + -) (C.l) 

7 
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Table CI. Parameters of the simulations. Here n is the number of layers in the 
Cayley tree, TV is the number of spins (the coordination number was fixed at z = 3), 
TVrun is the number of runs averaged over to extract error-bars, TVgwocp is the number of 
Metropolis Monte Carlo sweeps, Tmin and Tmax are the lowest and highest temperatures 
simulated, and TV^ is the number of temperatures used for parallel-tempering. 



and argued that as n — )■ oo, this system undergoes a sharp change in behavior at the 
critical temperature given by 7tanh(/3J) = 1. Now we notice that in fact, this quantity 
can show the transition even when n is finite. The key point is that at precisely the 
transition temperature, the quantity (7tanh(/9J))" becomes independent of n, therefore 
by plotting this quantity as a function of temperature for various system-sizes on the 
same graph, we can instantly obtain the transition temperature as the point at which 
the various curves take the same value. This is done in the main paper. Next we showed 
in the main paper that the quantity ^'^^ where Vt = -^/((croTVfn)^) — (co^n)^ is the 
RMS deviation, has a sharp change at the spin-glass transition temperature given by 
X = 7tanh {(3 J) = 1 in the thermodynamic limit, n -^ oo. Now, we point out that the 
mean square deviation Q"^ itself can be used in a finite-size-scaling study to obtain Tsg- 
We showed in the paper, based on the techniques described in the last section above, 
that, 

"1 77,-1-1 "1 

Q' = ((aoM„)2) - (aoMn)' = TV„(1 + ^^^ ^ - -x^). (C.2) 



7 X — i 7 



For finite n, it is better to write this as 



fi2 = ]v„(l + ^ — -x{l + x + x^ + --- x"-i) - -x"). (C.3) 



7 7 



It is clear that when x = 1, fi^ = Nnu"^^^. Therefore, the quantity -^^ becomes 
independent of system-size at this particular value of the temperature; thus plotting 
this quantity as a function of temperature for various sizes on the same graph, we are 
able to obtain the spin-glass transition temperature Tsg-, as shown in the main paper. 
Our simulation used the standard Metropolis algorithm in conjunction with parallel- 
tempering which helps speed up the equilibration. Table ICll collects the parameters of 
the simulations. Equilibration was verified by keeping track of the energy; those points 
(at very low-temperature, and large system sizes) for which the energy did not stay the 
same within error-bars over the last two sweeps, were discarded. 
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Table Dl. Parameters of the simulations on ± J spin glass. Here n is the number 
of layers in the Cayley tree, N is the number of spins (the coordination number was 
fixed at z = 3), A'samp is the number of samples averaged over to extract error-bars, 
TVswccp is the number of Metropolis Monte Carlo sweeps, Tmin and Tmax are the lowest 
and highest temperatures simulated, and Nt is the number of temperatures used for 
parallel-tempering. 



Appendix D. The Spin-Glass Transition Temperature 

Here we consider the ±J spin glass on the Cayley tree given by the Hamiltonian 

where Jij are independent identically distributed random variables that can take values 
± J with equal probability, and where the sum involves pairs of spins which are adjacent 
on the tree. Instead of the spin variables, like before, we choose the spin at the root of 
the tree o"o, and the bonds 9ij = a^Cj to define the configuration of our system. The 
6 variables can take values ±1, which make them effectively spin variables too. The 
Hamiltonian then takes the following simple form 

The partition function for this model is identical to that for the Ising ferromagnet since 
it is an even function in J, but more importantly, with the help of the same trick of 
inserting squares of all the spins that lie between any two given spins, we can now write 
down the correlator: 

{a,aj) = ±tanh{/3JY^^ = ±a'^^^ . (D.3) 

Whether a plus sign appears or a minus sign appears depends on the number of bonds 
between i and j that are ferromagnetic. We will now compute the 'Edwards- Anderson' 
type order-parameter with the knowledge that memory is the order-parameter for the 
ferromagnet, i.e we are interested in 

q = [{aoM^)X„g, (D.4) 

where [■ ■ ■]avg means an average over disorder. 
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Figure Dl. Monte Carlo simulation of the ±J spin glass no a Cayley trees with 
coordination number z = 3. The figure obtains the transition temperature popularly 
called Spin-glass Tsa HH] ■ While this temperature is of importance also in the regular 
ferromagnet as shown in the main paper, the reason for its name comes from studying 
the spin-glass as shown here. 



Computing the correlation between the root spin and the sum of the spins at the 



n level we obtain: 



(aoM„)= Y. K^.) = tanh(/3Jr J^ ±1, 



(D.5) 



Nnterms 



where Nn is equal to the number of spins in the n^^ layer: A^^ = z{z — 1)*^ ^, where z is 
the coordination number for the type of Cayley tree shown in Fig{Tl Squaring we have: 

(aoM„)2 = tanh(/3J)2"( ^ ±1)( ^ ±1) (D.6) 

Nnterms Nnterms 

The ±1 terms in the two sums are completely random, so when we square and average 
over disorder, only the 'diagonal' terms which give (±1)^ = 1 stay, all other terms 
vanishing because -|-1 is as likely as —1. Therefore, defining 7 = 2 — 1, we have 

1, 



[(aoM„)2],,g = (7tanh(/3J)2f(l 



7 



(D.7) 



It is clear that the behaviour of [(croM„)^] for 'long times' or, equivalently, in the 
'thermodynamic limit' n — )■ 00, changes completely depending on whether 7tanh(/3J)^ 
is > or < 1. So, 7tanh(/3J)^ = 1 signals a transition which is understood to be 
the so-called spin-glass critical temperature Tsg = -. — t-tttt-t- Fig. IDll shows that 



tanh-i{-^)- 



this phase transition can be obtained from Monte Carlo simulations of finite-systems. 
Again we exploit the convenient fact that at Tsg, the quantity (7tanh(/3J)^)" becomes 
independent of n, therefore by plotting g as a function of temperature for various system- 
sizes on the same graph, we can instantly obtain the transition temperature as the point 
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at which the various curves take the same value. Our simulation once again used the 
standard Metropolis algorithm in conjunction with parallel-tempering which helps speed 
up the equilibration. Table IDll collects the parameters of the simulations. We make 
the observation here that the phase below Tsa is not really a 'spin-glass' in the sense 
that there is no complex free-energy landscape. Two essential characteristics of a spin 
glass are frustration and disorder; here, we have only disorder, no frustration. Thus the 
transition Tsg is more like a regular ferromagnetic phase transition. 
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